Can electrostatic field lift spin degeneracy? 
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There are two well known mechanisms which lead to lifting of energy spin degeneracy of single 
electron systems - magnetic field and spin-orbit coupling. We investigate the possibility for existence 
of a third mechanism in which electrostatic field can lead to lifting of spin-degeneracy directly 
without the mediation of spin-orbit coupling. A novel argument is provided for the need of spin- 
orbit coupling different from the usual relativistic considerations. It is shown that due to preserved 
translational invariance spin splitting purely by electrostatic field is not possible for Bloch states. 
A possible lifting of spin degeneracy by electrostatic field characterized by broken inversion and 
translational invariance is considered. 



I. INTRODUCTION 

The field of spintronics"'^, a term coined by Wolf^, has 
attracted considerable attention in the past 20 years. 
This interest was initially triggered by the discovery of 
the effect of giant magnetoresistance (GMR)'^''* which 
now has established commercial applications. A paral- 
lel interest in semiconductor spintronics has arisen due 
to the proposal for a ballistic spin-transistor^ which has 
recently been realized^. A generalization of the device 
operating in diffusive regime has been proposed^ as well 
as other devices^. They depend on the spin splitting 
caused by the spin-orbit coupling through several differ- 
ent mechanisms usually classified in two big groups. His- 
torically the first group often referred to as Dresselhaus 
spin-orbit coupling^^^'^, or alternatively, bulk inversion 
asymmetry^^"^'^ (BIA), has as a necessary condition the 
broken space inversion invariance of a bulk crystal. It 



also appears in a modified form 



1,11,12,14,15 



in quasi two- 



dimensional (Q2D) structures grown in crystals lacking 
inversion symmetry. The second group referred to as 
either Rashba spin-orbit coupling^^^^® or structure in- 
version asymmetry^di,i2 (SIA) leads to spin splitting in 
systems lacking macroscopic inversion symmetry. In all 
these cases the inversion asymmetry is merely a neces- 
sary condition and not a sufficient one. The other neces- 
sary condition is taking spin-orbit coupling into account 
in various ways but usually through the folding down 
procedure^^"^'^'^^ from 8 band or 14 band models. Here 
we consider whether the breaking of inversion invariance 
of electrostatic field in several classes of systems can also 
be a sufficient condition for lifting of spin degeneracy. 

We examine a symmetry argument^^^^^'^^ due to 
Kittel"'^^. We provide novel derivation of the symmetry 
argument^^^^'^'^^ relying only on the commutation prop- 
erties of the time-reversal K and space inversion / op- 
erators with the translation Tr^ and spin operators Su. 
This is unlike the usual derivation^^ relying on explicit 
form of Bloch states. We show a novel argument for the 
need for introduction of a spin-orbit coupling term simi- 
lar to the argument motivating Maxwell to introduce the 
displacement current term in the Maxwell's equations. 
We note that the symmetry argument^^^^'^'^^ involving 



time-reversal and space inversion invariance does not re- 
quire the presence of a spin-orbit coupling term in the 
Hamiltonian but merely to take the spin degree of free- 
dom of the electron into account. The application of this 
symmetry argument, to a nonrelativistic model of elec- 
tron with spin moving in electrostatic field characterized 
by discrete translational invariance, naturally suggests 
the hypothesis that electrostatic field with broken space 
inversion symmetry can lead to lifting of spin degeneracy. 
We explore the hypothesis by perturbation method treat- 
ment and show that if discrete translational invariance is 
preserved electrostatic field can not lift spin-degeneracy. 
If both translational and space inversion symmetry of a 
perturbing electrostatic field are broken a possible con- 
tribution to spin splitting by electrostatic field is indi- 
cated. Physically this would be naturally explained since 
we know^^' that every non-accidental degeneracy stems 
from symmetries of the underlying system and in general 
reduction of symmetries leads to lifting of degeneracies. 



In Sec. II we introduce the theoretical models. The 
basic symmetry argument is presented in Sec. HI. The 
transformation properties of Bloch states in a model ne- 
glecting spin-orbit coupling but taking into account the 
spin degree of freedom are given in Sec. HI A. They are 
used in Sec. HI B to show the appearance of at least four- 
fold degeneracy in the spectrum as a consequence of the 
combination of time-reversal and space inversion invari- 
ance and in Sec. HI C to show the lifting of spin degen- 
eracy of Bloch states as a consequence of broken space 
inversion invariance. How the presented symmetry ar- 
gument differs from the usual one^^"^^'^^ is discussed in 
Sec. HID. In Sec. IV a new argument for the introduc- 
tion of spin-orbit coupling is presented different from the 
usual relativistic arguments^^'^^'^'^'^'*. We then explore 
the hypothesis formulated in Sec. HID by perturbation 
method treatment in Sec. V. In Sec. VI we discuss cer- 
tain aspects of the utilized model and methods and we 
present our conclusions in Sec. VII. 



II. THEORETICAL MODELS 



satisfies the Bloch theorem^^, 



We first consider two basic models of a nonrelativistic 
electron moving in a pure electrostatic field characterized 
by discrete translational invariance. We focus our atten- 
tion on ID models and 3D models of the triclinic crystal 
system. The two classes of models differ from each other 
by the properties of the electrostatic potential with re- 
spect to space inversion symmetry. In both models we 
take spin degree of freedom into account but neglect the 
spin-orbit coupling term. The eigenstates of both mod- 
els are Bloch states^^, which can be characterized by two 
quantum numbers: (i) crystal wavevector and (ii) spin in- 
dex which is the eigenvalue of the it-component S'„ of the 
spin vector operator S. Furthermore we employ the stan- 
dard Born-von Karman periodic boundary conditions. 

The first class of models represents ID systems or 3D 
crystals of the triclinic pinacoidal symmetry class. Its 
Hamiltonian takes the form 
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(1) 



This represents a nonrelativistic electron moving in a 
pure electrostatic potential 4>o{r) with potential energy 
Vo = — e(/)o(r) characterized by translational invariance, 
[Tr„, '/'o(f )] = 0, time-reversal invariance [K, (f>a{r)] — 0, 
and space inversion invariance [/, (/)o(t")] — 0. 

The second class of models represents ID systems or 
3D crystals of the triclinic pedial symmetry class in which 
the overall potential does not possess space inversion in- 
variance. A general potential not possessing space inver- 
sion invariance can be split into parts symmetric with 
respect to space inversion, that is one that commutes 
with the space inversion operator /, and a part antisym- 
metric with respect to space inversion, that is one that 
anticommutes with /. We consider the electrostatic field 
of the second class of models as made of the symmet- 
ric part (f>oif) and an antisymmetric one (/'(r), satisfying 
I4){r)I^ = —(j){r). The Hamiltonian takes the form 
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(2) 



While the Hamiltonian in such systems possesses discrete 
translational invariance [H,Tji^] — 0, and time-reversal 
invariance, [H, K] — 0, it is no longer invariant with 
respect to space inversion, [H, I] ^ 0. 

In order to emphasize that we have taken spin degree of 
freedom into account we have written the Hamiltonians 
(1) and (2) with the 2x2 identity matrix ctq. 

III. SYMMETRY ARGUMENT 

A. Transformation of Bloch States 

We consider a general Bloch state \k,Su) the properties 
of which are identical in the two models considered. It 



where we use the active convention^"'^^ for the space 
translation operator Tr^ . The time- reversal transformed 
state \k',s'^) ~ K\k,Su) of a Bloch state \k,Su) is 
still a Bloch state because the time-reversal operator 
commutcs^^ with the spatial translation operators Tr^ . 
Applying the time-reversal operator K to the Bloch the- 
orem (3), taking into account that [Tji^,K] = and that 
K as antilinear operator does not commute with com- 
plex scalars c but satisfies the idcntity^^ Kc = c*K, 
one obtains the identity rR^|A;',s^) — e'^'^'^"'\k',s'^). 
Comparing it with the Bloch theorem, TR^\k',s'^) — 
e^^'^ '^^\k' ,s'^), shows that k' = —k. Applying K to 
the relation S'„|fc,s„) — s„|A;,s„) and using the rela- 
tion KSu = —SuK, which follows directly from the 
definition^^ of the time-reversal operator K, one obtains 
SuK\k, s-u) ~ —SuK\k, Su)- Thus the X-transformed 
state \k',s'^) = K\k,Su) is an eigenstate of Su with 
an eigenvalue — Sm, therefore s'^ ~ —Su- Summarizing, 
the time-reversal operator K transforms a Bloch state 
|fc,s„) representing an electron moving with a wavevec- 
tor k and a spin pointing "up" the axis u into the Bloch 
state I — fc, — s,.), 



A|/c,s„y — I fe, Su) 



(4) 



representing an electron moving in the opposite direction 
with a crystal wavevector k and a spin pointing "down" 
the axis u. 

The space translation and space inversion operators do 
not commute but satisfy^° the identity 
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Applying the space inversion operator / to the Bloch the- 
orem, (3), using Eq.(5) and the fact that / docs not act 
on the phase factor e~^'''^" one obtains 



r-fl„/|fc,s„)=e-*^''"/|fe,s„). 



(6) 



Since the /-transformed Bloch state I\k,Su) satisfies the 
Bloch theorem it is still a Bloch state, but in general with 
different quantum numbers |fc',s^) = /|fc,s„). Com- 
paring Eq.(6) with the Bloch theorem T-ij^ |fc', s^) = 
gtk fli.|fc' s^^ one obtains that k' = —k. By definition^^ 
the space inversion operator I commutes with any com- 
ponent of the spin vector operator s. As a consequence, 
using the usual procedure applied above to Su\k,Su) = 
Su\k,Su) gives that s^ = Su- Therefore, the space in- 
version operator / maps a Bloch state \k,Su), describ- 
ing electron motion with crystal wavevector k and spin 
pointing in the direction of axis u, into the state 



/|fc,s„) 



rv. Si, 



(7) 



representing an electron motion with the same orienta- 
tion of spin but moving in the opposite direction with a 
wavevector k. Using the definition of the conjugation 
operator^^ C = KI and Eqs. (4) and (7), one obtains the 
action of the C on Bloch states, 



C|fc,s„) 



Ifc, 



.) 



(8) 



B. Spin Degeneracy 



Supposing the spectrum problem of the Hamiltonian 
Hf) with space-inversion invariant electrostatic potential 
(po (r) solved its eigenvalue-eigenvector problem takes the 
form of the identity 



Hn\k,Su) =El\k,Su) 



(9) 



The eigenvalues E^ ^ of Hq are labeled with quantum 
numbers fe, Su and the eigenstates of Hq and |fe, Su) sat- 
isfying the Bloch theorem possess all the properties of 
Bloch states, in particular Eq. (4) and Eq. (6). 

Applying the time- reversal operator K to Eq.(9) and 
using the time-reversal invariance of the Hamiltonian 
[iJ, K] = 0, we obtain 



HoK\k,s^) = El^^K\k 



(10) 



Due to the time-reversal invariance of H^ the two linearly 
independent states \k,Su) and | — k,~Su) = K\k^Su) 
correspond to the same eigenvalue of Hq^ which we 
now denote as E^^ = E^l ^^ = -^-fc-s„- The pairs 
of linearly independent states (|fc, s^), | — fc, — s„)) and 
(I — fc, Su), \k, —Su)) span respectively the 2D subspaces 
e'g and e'^g of the Hilbert space of the single-particle 
system. The 2-fold degeneracy of the cigcnenergies Eg 



and E°, 



E° 






= E^ _g , to which the subspaces 



eg and e'^g correspond, is a consequence of the time- 
reversal invariance of Hq and is the realization of Kramers 
degeneracy^-^'^^ in the system described by Eq. (9). By 
their construction the subspaces Eg and e'^g are invari- 
ant with respect to the time reversal symbolically writ- 
ten as Ke'^g = £j_g. Comparing the spectrum equation 

Ho\-k, -Su) = E'ii^_^J-k,-Su) for a state |-fc,-s„) 
with Eq. (10) allows us to express the Kramers degener- 
acy in the studied system in the form 



EL 



E-k.-s 



(11) 



Acting on the left of Eq. (9) with the space-inversion 
operator /, using the hypothesis [H, /] = and the result 
| — fc, Su) = I\k, Su) from Eq. (7), one obtains the identity 



Ho\ -k,Su) ^ EL^J -k,Su) 



(12) 



It shows that the space-inversion invariance of the Hamil- 
tonian Hq requires that the Bloch states |fc,s„) and 
I — fc, Su) belong to the same eigenenergy. 



Ek = E^ 



E\ 



(13) 



Expression (12) shows that the states |fc,Su) and 
I — fc, s„), respectively belonging to the subspaces eg and 
e_e, rnust belong to the same degenerate eigenvalue E^ 
as a consequence of the space inversion invariance of iJg • 
Taking into account the consequences of the time-reversal 
invariance of Hq given in Eqs. (10) and (11), all eigenen- 
ergies E^ of Hq must be 4-fold degenerate. To every en- 
ergy value E''1 corresponds a four-dimensional subspace 
e°, which is a direct sum, e° = e^ + e^_g, of the two sub- 
spaces Eg and s'^_g. The subspace e^ is invariant and re- 
ducible with respect to space inversion and time-reversal 
written symbolically as Ke^ — e^ and Ie\ = e^. 

Combination of the space inversion and time-reversal 
invariance of the Hamiltonian Hq of the considered trans- 
lational invariant system is equivalent to spin degener- 
acy. Formally this is illustrated using the conjugation 
operator C = KI which commutes with the Hamiltonian 
Hq, [Hq^C] = if it commutes separately with K and 
/. Using the usual procedure of applying the operator 
C to Eq.(9) and taking into account Eq.(8) one obtains 
the identity Ho\k,—Su) = E'^ ^ |fc,— s„). Therefore the 
Bloch states |fc,s„) and |fc,— Su) describing an electron 
with opposite spins belong to the same degenerate en- 
ergy value E^. The subspace e° corresponding to E^ is 
invariant with respect to C . The spin degeneracy can be 
viewed also as a consequence of the SU{2) invariance of 
the Hamiltonian Hq. 



C. Broken Spin Degeneracy 

The problem for the spectrum of the Hamiltonian H 
shown in Eq (2) is given by the identity 



7/1/4,0-,, 



Ek 



(14) 



The Bloch states \k,(Ju) are the common set of eigen- 
states of the commuting operators H and Tr^ and E^u^ 
arc the corresponding eigenvalues of the Hamiltonian H . 
They transform among each other according to relations 
(4), (7) and (8) because the Hamiltonian H given in 
Eq (2) is translational invariant. 

The spectrum (14) of the Hamiltonian H possesses the 
time-reversal induced properties derived from Eq. (10) 
and Eq. (11). It is Kramers degenerate and each of its 
cigcnenergies are two-fold degenerate, i?K,,cr„ — £'-k,-<t„- 
However because the space inversion invariance of the 
electrostatic potential (t>{r), and therefore of H is bro- 
ken, the degeneracy due to space inversion invariance is 
lifted, i?K,(T„ 7^ -E_K,CT„ • This requires the hfting of spin 
degeneracy 



Ek 



^ Ek,,-^^ 



(15) 



The detailed proof follows in the next paragraph. 

The broken spin degeneracy is proved by applying the 
conjugation operator C to Eq. (14). However, because 



[H, /] ^ and hence [H, C] ^ wc can not interchange 
the positions of C and H. Instead, using 1 — C~^C and 
I**, —(Tu) = C\K,au) one obtains from Eq. (8) that 



CHC ^\K,-a-u.) ^ E^^^jK,-a^ 



(16) 



The Bloch state |a«, — cr„), as an eigenfunction of Tr^, 
and because [H,Tfi^] = 0, is stiU an eigenstate of H 
with an eigenvalue E^^^u^. However, because of broken 
inversion invariance, 



CHC- 



ikHk~^r^ = inr^ ^ h 



\i^,^(^u) is not anymore an eigenstate of H with the 
eigenvalue Ei^^^r^. Instead the Bloch state |k, — CTu) is 
eigenstate of some other operator H' — IHI^^ with 
the eigenvalue £^K,,cr„- Therefore the states |k, CTu) and 
\i^, —(Ju) do not correspond to the same energy. Elec- 
trons in Bloch states characterized by the same wavevec- 
tor K but having opposite spin orientations do not posses 
the same energy, the result given in Eq. (15). 



D. Discussion of Symmetry Argument 

In the well known treatment^'^^'^^'^^ it is supposed that 
the spin-orbit coupling is part of the model Hamiltonian 
and that the spin degeneracy is lifted by the spin-orbit 
coupling term 



Hso — 



Am^c^ 



& ■ {\7V{r) X p) 



(17) 



if the electrostatic potential V{r) — Vo{r) — e(t){r) does 
not possess space inversion invariance. However, close 
examination of the symmetry analysis developed in the 
text above shows that there is no such requirement. The 
symmetry analysis is valid also for a nonrelativistic model 
that contains just electrostatic fields (j)o{r) and (j){r) and 
does not contain spin-orbit coupling or magnetic field. It 
suggests that electrostatic field alone without the media- 
tion of spin-orbit coupling can lead to lift of spin degen- 
eracy given that (t>o{r) and (j){r) possess discrete transla- 
tional invariance and are characterized by preserved and 
broken space inversion symmetry respectively. 



IV. NOVEL ARGUMENT FOR 

INTRODUCTION OF SPIN-ORBIT COUPLING 

TERM 

Taking spin into account suggests examining for SU(2) 
symmetry. Since we consider models of a nonrelativistic 
electron in electrostatic field neglecting spin-orbit cou- 
pling, the Hamiltonians Hq and H commute with ev- 
ery component of the spin vector operator S and there- 
fore the Hamiltonians Hq and H are SU(2) invariant. 
They commute, [R'^{4>),H] = 0, with every spin rotation 
operator^^ Rui4>) = e"*"^*^" for arbitrary axis u and angle 



of rotation (/). As a consequence of the SU(2) invariance, 
the spin degeneracy for the model with Hamiltonians Hq 
and H must be preserved, in particular for H 



Ek 



Ek.-h 



(18) 



This result, however, contradicts the result (15) stem- 
ming from symmetry analysis based solely on time- 
reversal and broken space inversion invariance when spin- 
orbit coupling is neglected. 

A possible way to resolve this inconsistency is the intro- 
duction of a term in the model Hamiltonian that breaks 
the SU{2) invariance and at the same time is consis- 
tent with the different cases of symmetry analysis in- 
volving just time-reversal and space-inversion operators. 
The spin-orbit coupling term (17), which docs not com- 
mute with any spin-rotation operator^^ _R|j(0), satisfies 
the above requirements and resolves the noted inconsis- 
tency. We interpret this as a novel argument for intro- 
duction of the spin-orbit coupling term different from the 
usual^'^'^'^'^'* purely relativistic considerations, which give 
an incorrect numerical factor by 1/2. This difference is 
accounted for by Thomas precession or by taking the 
nonrelativistic limit of the Dirac equation^^'^^'^'*. Thus a 
realistic model of electron dynamics requires taking into 
account the spin-orbit coupling term as a minimum; oth- 
erwise we would encounter the above mentioned incon- 
sistencies. Therefore in all subsequent models treated 
within perturbation method the spin-orbit coupling term 
is part of the considered Hamiltonian. As a consequence 
SU(2) symmetry is always broken and there is no require- 
ment for preservation of spin-degeneracy. 



V. PERTURBATION METHOD TREATMENT 

A model with spin-orbit coupling does not exclude the 
possibility that electrostatic field by itself leads to lift- 
ing of spin degeneracy suggested by the symmetry anal- 
ysis based on time-reversal and space-inversion invari- 
ance in Sec. HID. It merely shows that if such splitting 
exists it will lead to additional numerical factor in the 
spin splitting already caused by spin-orbit coupling when 
the space-inversion invariance of the electrostatic field is 
broken. Since the crystal wavevector k varies in discrete 
steps we investigate this option using standard stationary 
perturbation method for degenerate levels. 

The unperturbed Hamiltonian is Ho = ^ + 
Voir), where Voir) possesses space inversion invariance, 
[Vb,/] = 0, while the perturbation SV consists of an 
electrostatic potential (j>{r) which is odd with respect to 
space inversion, {0(r),/} — 0, and the spin-orbit cou- 
pling terms 



f/2 
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;& ■ (VVoir) X p) , (19a) 

a ■ (V</.(r) X p) . (19b) 



Unlike in the standard k ■ p method where the fe = 
stationary states are used as unperturbed basis, we use 
the stationary states \k,s-u) of Hq for arbitrary k ^ 0. 
This choice is naturally suggested by the symmetry anal- 
ysis above since for A; = we have just two-fold Kramers 
degeneracy that is not lifted as far as time-reversal invari- 
ance is preserved. We consider ID model or 3D model of 
triclinic pedial and triclinic pinacoidal systems in which 
cases the degeneracy of every energy level E'^ of the 
unperturbed Hamiltonian Hq is exactly four-fold. The 
corresponding subspace e^ is spanned by the four Bloch 
states |A;, Su), \k, — «„), | — fe, s„) and | — k, — «„). 

The first-order correction Ej^ to the energy eigenvalue 

Ej, and the zeroth-order states |0) arc determined from 
the eigenvalue equation 



p,"<5yp°|0fe) = ii;W|0fc), 



(20) 



where P° is the projector to the subspace e^ and 6V — 
—e(j){r) + Ui + f72.The first-order correction EJ, ' to 
the energy is determined by the solution of the secular 



equation det 
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E 
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corresponding to the 



eigenvalue problem (20). The condition for time-reversal 
invariance simplifies the secular equation by introducing 
relationships between its matrix elements, 



-E^ 
c* 
d* 




(1) 
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6-K 



(1) 




d 
— c 

-E, 



(1) 



= , (21) 



where the matrix elements a = 01+02+0:, b = 61+62+0, 
c = C1+C2 and d = di+d2+P have contributions from the 
three different perturbing terms. The matrix elements 
due to the spin-orbit coupling term Ui are denoted as 
ai, 61, ci, di, the ones due to U2 as 02, 62, C2, d2, and 
the ones due to the electrostatic field <j){r) as a and /?, 

oi = (s„,fe|L/i|fe,s„) , 02 = (s„,fe|f72|fe,s„) , (22a) 
61 = (-s„,fe|L/i|fe,-s„) ,62 = (-s„,fc|C/2|fe,-s42?b) 
ci =: (s„,fe|f7i|fe,-s„) , C2 = (s„,fe|f72|fe,-s„) (22c) 

di = {su,k\Ui\ -k,Su) , ^2 = {su,k\U2\ -fe,s„;(?2d) 
a = {k\SV{r)\k) , (3 ^ {k\5V{r)\ - k) . (22e) 

Taking into account the properties of the electrostatic 
potential (p{r) and the terms Ui and U2 leads to further 
simplification of the above result. The inversion invari- 
ance lUiI^ ~ Ui requires that oi — 61, ci = and that 
di is a real number. The condition that the term U2 
must be odd 11/2!^ = —U2 with respect to space inver- 
sion invariance, requires 02 = —62 and d2 to be purely 
imaginary, while it does not place any restrictions on 02- 
The fact that the electrostatic field (/)(r) is odd with re- 
spect to inversion invariance, I(j>{r)I~^ ~ — (/)(r), requires 
a = —a and thus a = 0, while /3 needs to be purely 



imaginary. Taking into account these simplifications the 
first order corrections to the energy take the form 



E 



(i)± 



— Ol ± 



C2 
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(3\' 



(23) 



Eq. (23) shows that the spin-orbit coupling term Ui 
possessing space inversion invariance leads to a shift 
of the first-order energy levels by the amount oi ~ 



A. Perturbation with Translational Invariance 

The expression (23) shows that the four-fold degen- 
eracy of E'^ can be lifted up to two twofold degenerate 
levels in first-order perturbation method. The spin split- 
ting caused by the spin-orbit coupling term U2, which 
lacks inversion invariance, is embodied in the matrix ele- 
ments 02, C2 and 62 within this approach. The presence of 
the matrix element /3 in the energy correction (23) leaves 
open the option that the electrostatic field leads to an ad- 
ditional contribution to the spin splitting. However, since 
the condition /(/)(r)/+ = —(j){r) requires that either (3 is 
purely imaginary or 0, the consideration of time- reversal 
invariance and broken space inversion invariance within 
perturbation method does not offer conclusive evidence 
that /3 ^ 0. 

We have supposed that the nonperturbed and per- 
turbed electrostatic potentials possess discrete transla- 
tional invariance, Tji^(f>{r)T^ — '/'(^)- Consequently, 

the terms Ui and U2 also possess discrete translational 
invariance, Tji^UiT^ = Ui, and Tji^U2T^ = U2 
respectively. Using this and the Bloch theorem for 

^2 = (su,fe|^2| - k,Su) we find 

(s„,fc|C/2| -fe,s„) - {s^,k\f+^U2fRj-k,s^} = 

^e^'''-''{s^,k\U2\-k,s^p4) 

which for arbitrary fe and translation vector H can be 
satisfied only if ^2 = 0. Similar reasoning shows that 
di ~ and /3 — 0. Thus the first-order energy correction 
(23) takes the form 
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(25) 



Therefore it is not possible to have lifting of spin degen- 
eracy of Bloch states purely by translationally invariant 
electrostatic field, which is odd with respect to space in- 
version invariance, but the reason for this is the discrete 
translational invariance of the perturbing electrostatic 
field (t>{r). This is interesting since in the usual sym- 
metry argument translational invariance is not explicitly 
considered, nor any attention is devoted to it. Further- 
more, all energy corrections of higher order disappear due 
to translational invariance since all (sit, fc|C/|fe', s„) — 0. 
This makes Eq. (25) exact to all orders in the perturba- 
tion method treatment. 



B. Perturbation without Translational Invariance 

Natural continuation of the above line of thought is to 
consider perturbing electrostatic field 4>'{r) with corre- 
sponding energy V'{r) = — e0'(r) not possessing a center 
of inversion, /V'(r)/+ 7^ Vif), and lacking translational 
invariance 



TR^V'ir)f+ + V\r) , 



(26) 



to which we will refer as external. Its corresponding first- 
order matrix element which is not necessarily zero due 
to time- reversal invariance is /3' = (sM,fc|T^'| —k,Su) = 
(— s„,fc|y'| — k,—Su). It is defined between zeroth- 
order states with spin, but since the electrostatic field 
does not contain operator acting directly on the spin- 
degree of freedom it can be written shorthanded as 
P' = {k\V' \ — k) . Because V'{r) is not invariant with re- 
spect to space inversion, the disappearance of the matrix 
element /?' is not guaranteed. Examples of electrostatic 
fields possessing these properties are the built-in poten- 
tial in p-n junctions and externally applied electric fields. 
These are usually orders of magnitude smaller than the 
bulk electrostatic fields. 

By using the Bloch theorem and the inequality (26) we 
obtain the inequality 



(s„,fe|T/'(r)|-fe,s„)^ 



2ikR, 



s„,fc|T/'(r)|-fc,.s„) 



(27) 
for the matrix element /3'. Therefore for any given 
k and arbitrary translation vector R we must have 
P' ^ e^'^^'^P' . This condition excludes the possibility the 
matrix element to be equal to zero and thus /3' 7^ 0. 

In calculating the first-order corrections to the en- 
ergy we neglect the matrix elements a', b' and 
c' stemming from the spin-orbit coupling term 
^J^i^-i O- ■ {VV'{r) X p) since they should be much smaller 

than /3' = {su,k\V'\ — k,Su) because of the prefactor 
.22 - Since /3' ^ it will appear in the expressions 
for the energy splitting. For the triclinic pedial system 
in which the bulk potential does not possess space in- 
version invariance, the first-order corrections to the en- 
ergy, Ef,' = ai ± y^a^ + |c2p -I- l/3'p, contain contribu- 
tions from the spin-orbit coupling 02 and C2 and a contri- 
bution purely from electrostatic field in /3'. The original 
fourfold degeneracy is lifted to two distinct two-fold de- 
generate levels. The energy splitting between them is 
given by 



a4'^ = 2 



|C2| 



\I3' 
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(28) 



For ID model or triclinic pinacoidal system where the 
bulk electrostatic potential possesses space inversion in- 
variance, a2 and C2 are identically zero. The expres- 
sion for the first-order energy corrections takes the form 



E 



(i)± 



ai ± 1/3' I . The energy splitting. 



between the two twofold degenerate levels has contribu- 
tions only directly from the electrostatic field through 



C. The lifted degeneracy 



The energy splittings (28) and (29) can be interpreted 
as pure spin-splitting for the zero-th order Bloch states. 
This is easily demonstrated using the invariances of the 
subspaces e^,.^ as shown in the following paragraphs. 

Since the fourfold degeneracy of Ej^ ' is not completely 

(1) 

k± 



removed, each of the first-order corrections to energy E^ 
corresponds^^ to a two-dimensional subspace e)^^ of the 
four-dimensional subspace e\. of the unperturbed prob- 
lem (9). The subspaces e^,.^ are mutually exclusive and 

(0) , JO) _ (0) 



their direct sum e^"^ -I- £^,_ 

-(0) 



£^ is the four-dimensional 



subspace e^, . The zeroth-ordcr state |0a;), which is the 
projection of the perturbed state |^) onto the subspace 

ejj, , cannot be determined uniquely, only its belonging 

to one of the subspaces e^^ can be inferred from Eq. (20). 
Symmetry arguments suggest that Kramers degener- 
acy is preserved and therefore the remaining twofold de- 
generacy in first order is the Kramers degeneracy. In or- 
der to show this wc suppose the equation for first-order 
energy correction to be solved and consider it as identity 



PfV'(r)Pr|0±).42lO^) 



(30) 



AE, 



(1) 



2|/3'| , 



where the states |0^) are arbitrarily chosen zeroth-order 
states belonging to the subspaces el.J., respectively, and 
therefore to ej. . Applying the time reversal opera- 
tor to Eq.(30) and using that [K^PJ^ '] = 0, because 
e° is invariant with respect to K, [K,V'{r)] — by 
hypothesis and Ef.J. is real, we obtain the identity 
P^°V'(r)pf 'X|0±) = e12k\0^). The ^-transformed 
state K\0^) is orthogonaP^ to |0^) since we consider a 
single-electron system taking into account the spin degree 
of freedom . As shown above, K\0^) also identically sat- 
isfies Eq. (30) with the eigenvalues £^^.1: and therefore be- 
long to the subspaces £^^7 respectively. Since we choose 
|0 ) arbitrarily, and by the above argument every K\0^) 
belongs to ej.]., the subspaces ej^J. are invariant with 

respect to the time-reversal operator K, Ke^.^ = ei.^. 
Thus the remaining two-fold degeneracy in first order is 
precisely the Kramers degeneracy. 

The two-dimensional subspaces £^2 are, however, not 
invariant with respect to the space- inversion operator be- 
cause the perturbation V'{r) does not commute with it, 
[/, V^'(t)] 7^ 0. This is proved by again applying the time- 
reversal operator / to the identity (30) satisfied by arbi- 
(29) trary state |0^) G £^°2- However, since [I,V'{r)] ^ 



and by using, I ^I = I, Eq.(30) transforms into the 
identity, pf ^/y'(r)/-ip^"'|0±) = ^f^^lO"^) for the /- 
transformed state /|0^). This shows that /|0^) does 
not satisfy Eq.(30) with eigenvalue E^J. but a different 
equation with the eigenvalue Ef^J. because V ^ IV' I^^. 
Therefore the /-transformed states /|0^) do not belong 
to the subspaces ej^J_. The choice of |0''=) is arbitrary 
apart from the condition |0^) G el.^ and therefore for 
every state |0^) belonging to e^_^ the /-transformed 
state /|0^) does not belong to sj^^. However ,the four- 



(0) 



-(0) , ^(0) 



-(0) 



e^"^ is invariant with 



dimensional subspace e 

respect to /, /e^, = e^, and therefore /|0^) must be- 



(0) 



long to e^ . Since /|0^) does not belong to sj^J. the 
only remaining option is that it belongs to the other 

two-dimensional subspace sl.^- Thus the space inver- 
sion operator / maps every state |0^) G e^^ to a state 
/|0+) £ el,_ belonging to the other two-dimensional 
subspace e^_ and vice versa, symbolically written as 



,(0) 



-(0) 



Is^^J_ = e)^'' In other words the zeroth-order states |0^) 
and /|0^) belong to the two different first-order energy 
corrections Ef,J_ separated from each other by the energy 
difference A/^j, '. 

Now consider any arbitrary zeroth-order state with 
some spin polarization belonging to the subspace e^,! 
which we denote as |0+} and apply the spin- flip 
operatori9 C to obtain C\0+) = ik\0+) = /|0+'). From 
the previous results we know that |0+') — /f |0+) also be- 
longs to the subspace e^ / , while the /-transformed state 
(710+) = /|0+') belongs to the other two-dimensional sub- 
space e5j,_. Therefore the spin- flip operator C, Eq. (8), 
maps any zeroth-order state |Ofe) e el.^ to a state be- 
longing to the other two-dimensional set \0k) £ e^^^ sim- 
ilarly to the space inversion operator /. By definition a 
zeroth-order state corresponds to the subspace e^.^ if it 
satisfies Eq.(20) with the corresponding eigenvalue Ej^J.. 

So the zeroth-order states |Ofe) G ej^J_ and C\Ok) £ £^,2 
characterized by identical quantum numbers but describ- 
ing opposite spin orientations belong to the two different 
first-order energy corrections Ef,J_. This constitutes the 
proof. 

A question might arise as to why should a matrix el- 
ement /3' which is diagonal in spin indices be thought of 
as a spin splitting even just for the zero-th order states. 
A perspective to understand the issue is that indeed if 
the secular problem to be solved reduces to 2 x 2 matrix 
then matrix element diagonal in spin indices can not be 
interpreted as spin splitting. However we consider a 4 x 4 
matrix and the matrix element /?' while diagonal in spin 
indices is not diagonal in the 4x4 matrix. They are on 



the third diagonal and they lead to lifting of the origi- 
nal four-fold degeneracy into two two-fold degeneracies. 
Above we have presented detailed argument for the in- 
terpretation of these as pure spin splitting for the zeroth 
order Bloch states. 

However the perturbed states |^) arc no longer Bloch 
states because the perturbation V'{r) breaks the space 
translation invariance. As a consequence the crystal 
wavevector k is no longer a good quantum number for 
the perturbed states and in higher orders in the perturba- 
tion method the lifted degeneracy can not be interpreted 
purely as spin degeneracy. Instead it can be interpreted 
as lifting of some sort of orbital-spin degeneracy similar 
to the fine structure splitting due to spin-orbit coupling in 
atomic systems. This is within the perturbation method 
treatment. 



VI. DISCUSSION 

The primary goal of the study has been to explore the 
hypothesis whether an electrostatic field can lift the spin 
degeneracy of Bloch states. Integral new contributions of 
this have been the novel point for introduction of spin- 
orbit coupling presented in Sec. IV and the null result 
of Sec. VA, which have naturally shaped the presenta- 
tion. The hypothesis presented in mathematical detail 
in Sec. HI and discussed in subsection HI D has been in- 
spired by a symmetry argument which to the best of our 
knowledge is due to Kittel"'^^. The original argument has 
been used to predict the lifting of spin degeneracy in the 
presence of spin-orbit coupling. In subsection HI D we 
noted that it is valid also for a model of electron moving 
in an electrostatic field in which we take spin degree into 
account but neglect spin-orbit coupling. Such a model is 
easily justified on the ground that in nature there is no 
electron without spin. 

The trivial consequence of such a model would be dou- 
ble degeneracy of all levels due to the spin-degree of free- 
dom. In such a model, where spin-orbit coupling is ne- 
glected, SU(2) symmetry is preserved as noted in Sec. IV. 
Thus the double spin degeneracy of such a model can be 
viewed also as a consequence of SU(2) symmetry. Of 
course when spin-orbit coupling is taken into account, 
SU(2) symmetry is broken^^, and it is not possible to 
introduce spin quantum number as quantum number of 
type constant of motion. However we treat spin-orbit 
coupling terms as perturbations and in the zeroth-order 
model the spin can be introduced as a separate quantum 
number. 

The original symmetry argument due to Kittel^^ does 
not use the irreducible representations of the point 
groups. We considered formulating the symmetry argu- 
ment using the irreducible representations of the point 
groups, but for the particular purpose we concluded that 
working with the symmetry operators themselves is suffi- 
cient. This is so, because we consider a conceptual ques- 
tion, therefore we test the hypothesis on the simplest pos- 



sible systems which posses the required characteristics of 
the problem. These are ID model and 3D crystals of the 
triclinic crystal system: the triclinic pedial and triclinic 
pinacoidal crystal classes. These are the least symmetric 
classes, which have in common only a rotation by 27r and 
a time-reversal symmetry. In addition the pedial class 
is characterized with a broken space inversion symmetry, 
while the only other symmetry element of the pinacoidal 
class is the space inversion. 

The compatibility relations of the irreducible repre- 
sentations can be used to calculate the effects, including 
lifting of degeneracy, by electric fields which do not break 
the translational invariance. However, this approach 
would not offer an answer to the hypothesis whether an 
electrostatic field alone would break the spin degeneracy. 
This is so because based solely on the symmetry argu- 
ment wc cannot say whether the spin splitting is caused 
by the electrostatic field alone or by the spin-orbit cou- 
pling. This is the reason to use the perturbation method. 
Indeed, when we use the perturbation method it turns 
out that if the translational invariance of the electrostatic 
field is preserved no spin-splitting occurs. On the other 
hand the perturbation method treatment indicates a lift- 
ing of the spin degeneracy when both the translation and 
inversion invariances are broken. However, if the trans- 
lational invariance is broken we can no longer talk about 
space symmetry groups and their subgroups - the point 
groups. This is so because the discussion of the space 



groups requires the translational invariance of the crys- 
tal lattice. 



VII. CONCLUSIONS 

We have scrutinized the symmetry argument based on 
space inversion and time-reversal invariance predicting 
the appearance of spin splitting in case of broken space 
inversion symmetry. A novel argument for the need for 
introduction of a term breaking SU(2) invariance like the 
spin-orbit coupling term has been presented. This ar- 
gument is different from the usual arguments based on 
special relativity used for the introduction of spin-orbit 
coupling term. We have shown that in systems possess- 
ing discrete translational invariance it is not possible to 
have spin splitting solely by electrostatic field with bro- 
ken space inversion symmetry due to the preserved trans- 
lational invariance while a spin splitting exists due to the 
combination of spin-orbit coupling and electrostatic field. 
The possibility for lifting of spin degeneracy due to elec- 
trostatic field without the mediation of spin-orbit cou- 
pling has been investigated using perturbation method 
treatment suggesting its possible existence in systems 
characterized by both broken space inversion invariance 
and translational invariance. 
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